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Abstract 



A method is established which allows the calculation of the walk dimensions for Sierpihski- 
type multifractals. The multifractal scaling behaviour of the average time needed to cover 
a distance in the Sierpihski-type multifractals is shown. We calculate the Renyi dimensions 
for the average-time-multifractal and apply the /(a)-formalism. 

1 About the objects we are considering 

The topic of our interest is the scaling behaviour of a certain set of geometrically constructed 
multifractals, the so-called Sierpihski-type multifractals. Strictly speaking, we generalize the set 
of Sierpihski-type fractals introducing a supplementary structure. 

First some conventions. Any pure Sierpihski fractal is characterized by the embedding Eucli- 
dian dimension d e and the base b. Scaling the space coordinates with the factor b yields a fractal 
point-set, which is indentical to the original fractal point-set. Further, we will call the simple 
geometric object representing the first step in the construction of a pure Sierpihski fractal an 
elementary generator of the base b. 

There are two ways to construct Sierpihski fractals. On the one hand, we can manage it 
starting from a <i e -dimensional simplex. We divide all edge lines in b parts. By connecting the 
division points with parallels to the edge lines we get a set of disjunct subcells, also simplices. 
Now we wipe out all simplices whose orientations differ from that of the original simplex. The 
number of remaining subcells in an elementary generator can be computed by the recursive 
formula: 



In the following we restrict ourselves to the cases d e > 2. In the above picture we approach 
the Sierpihski fractal from the nearest higher-dimensional simplex. But before finishing the 
construction we are not able to define a reasonable model of walk. On the other hand, we can 
compose it as a gasket of one-dimensional line pieces. For every construction step we can define 
a walk model on the gasket. After an infinite number of repetitions both constructions yield the 
same fractal. But for our intentions the second picture will be better. 
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or their explicit equivalent: 





Figure 1: After the second construction step; here we have 2 generator types (bases 2 and 3) 
with the ratio of average numbers 2:1 

Our Sierpihski-type fractals result by mixing some elementary generators of different bases 
(figure 1). After numbering all elementary generators we define a matrix A, whose elements 
ciij are the average number of i-type elementary generators in j-type elementary generators. 
Sometimes we will use the factorized form: 

dij = gijiij with rij = n(bj,d e ). (3) 

Here §ij is the probability to find subcells with i-type elementary generators in the j-type parent 
elementary generator. Obviously we have 

$>j = l,Vj. (4) 

i 

The knowledge of the matrix A does not completely determine our multifractal's geometrical 
structure. Also there is a stochastic element in it; we can not say anything about the exact 
distribution of subcells including different elementary generators in any parent elementary ge- 
nerator. The only thing we are sure is, that after a finite number N of construction steps 
we can realize the distribution coefficients gij as exactly as we wish. Thus, we will neglect in 
our calculations all contributions of the first N construction steps and treat the multifractal as 
completely characterized by the given matrix A. 

2 Renyi dimensions and /(a)-formalism 

To compute the Renyi dimensions D q and to apply the /(a)-formalism we ought to mention 
some facts about the notion of dimensions for multifractal point-sets (see [1]). The so called 
local dimension a(r) describes the local scaling behaviour of some extensive (physical) property 
p of our multifractal system. Marking all the points with the same value of a yields another 
fractal point-set. Its Hausdorff dimension is f(a). The Renyi dimensions are a set of dimensions 
parameterized by q. This set includes the Hausdorff dimension Dq, the information dimension 
D\ and the correlation dimension D<i- These dimensions D q are so called global dimensions. 
That is, they characterize the whole multifractal point-set and do not have any local notion. 
But they are closely related to a and f(a) by the following formulas: 
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and: 

D q = ^ 1 (qa(q) -/(«(«))). (6) 

To get the Renyi dimensions of a given multifractal it was defined a partition function for any 
splitting {S{} of our multifractal into K cells: 

r(,,r,{5,},0=Ef ■ (7) 

i=l i 

Here pi is the value of the considered property and /; is the edge-length of the i-th cell. In the 
limit of infinite small cells the sum above converges to T(q,r) for any sequence of {Si}. It was 
been shown that demanding T(q, r) to be of order 1 (exactly speaking: < T(q,r) < oo) leads 
to 

r(q) = (q-l)D q . (8) 

Thus we get 

f(a(q)) = qa(q)-r(q). (9) 

Now we will return to our multifractals. Here we use the set of our fractal's own subcells in the 
actual construction step as {Si}. Just the question arises: What happens with the above sum 
when continuing the construction by one step? Obviously we must treat all contributions to the 
above sum from subcells including the same elementary generator (after the next construction 
step) in the same way. Therefore it seems useful splitting the sum T(q, r) into the elements of a 
vector 7(</,t). Each element of this vector is directly related to one elementary generator type 
incoming in our multifractal's construction; it contains the contribution of all subcells including 
the same elementary generator type in the next construction step. 

It is necessary to guarantee the right treatment of these partial sums in the further construc- 
tion. It is the moment to use the matrix A. In each construction step it sorts the next sum in 
the correct way and places all terms in the right partial sum. The change of the summands due 
to the scaling will be managed by the two matrices pij and kj. pij is the amount of the extensive 
(physical) property of the subcell i placed into the cell j with 

9ijnjPij = 1 , Vj . (10) 

i 

The lij = lj = 1/bj are the edge-lengths of the i-type subcells in the parent cell of the type 
j with the edge-length 1. Combining all this we get the transformation rule of j(q, r) for the 
transition to the next construction step: 

^+i) = M1 w with M = (m .. )= L.zL] . (ii) 

The matrix M contains all needed information to compute the multifractal dimensions. It 
only misses out the exact geometrical structure. That is, we know nothing about the concrete 
distribution of a(r). However, in M we find all information about the global density of cells (in 
the limit: points) with the same value of a and that is sufficient for our calculations. Therefore 
it seems reasonable to name M the "multifractal-generator(-matrix)". In our description "to 
be of order 1" is transferred from T to the largest eigenvalue of M. This is equivalent to the 
demand for the existence of a nontrivial matrix M°° with non-negative elements. To make sure 
that its largest eigenvalue must be equal to 1. It can be proven easily. 



3 



The walk dimension d 



to 



For the following we need the time scaling factors A 1 for each type of elementary generators 
incoming in our multifractal's construction. Thus we start with proposing a method for their 
calculation. Its explicit derivation shall be omitted here. First we mark one of the corner points 
as the starting-point and the d e other corner points as end-points for any way of particles passing 
through our elementary generator's network. In our simple walking model the time needed to 
cover the distance between two nearest neighbor points of the elementary generator has the 
value 1 and there is no waiting time between two walking steps. Now we number all elementary 
generator points giving the same number to all that points in equivalent positions relating to 
the starting-point. Further, all end-points shall absorb the arriving particles completely. The 
next step is to write down the propagator containing all transition probabilities between nearest 
neighbor points. This suggests to use a vector-matrix-notation. The i-ih element of the vector 
X^ used here contains the probability to find the particle located at the i-th point after the k- 
th step. The application of the propagator matrix Q = (qij) creates the probability distribution 
after the next step: 

${k) = Q k$(f>) (12) 

The contribution to Q from the end-points can be neglected due to their absorbing property. 
This makes Q a little smaller. Further, it is reasonable to define the vector X: 



k=o \k=0 

The absorbing end-points yield 



£qM*(°). (13) 
\k=0 ) 

Q°° = 0. (14) 



Therefore we are allowed to write 



E = = " G)" 1 ■ (15) 



Here E is the identity matrix. We get 

-? = ^-?<°>. (16) 

If h{ is the number of equivalent points specified by the same index i, the probability of passages 
through the link ij in the k-th step is 

,,s ( h;x^ + h:x^\;: for i ^ j 

EM(i~j)={ ^ ' > *' TJ (17) 

{ hix) 'qji for i = j . 

This yields the total (average!) number of passages through the link ij counted over all possible 
ways: 

I hiXiqa + hiXiqa for i ^ j 1 

E(i~j) = \ ; ^ J jqij t T . (is) 

[ hiXiqji for i = j . 

Now we get 

A = - S(i ~ j) = lij^jhj . (19) 



1 For any diffusion (a; 2 ) (xt d ™ and d w = 



In \ 
In b 
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For having neglected the absorbing end-points as well in X^ and X as in the propagator matrix 
Q, we are allowed to simplify the above formula to 




iXi = H T X = H 



1 



(20) 



E-Q 



Back to our principal problem the question arises: Can we use the method described in the 
preceding section to compute the walk dimension d w ? Yes, we do as will be seen soon. 

The following comments should enlighten our idea. In our model the walking particles leave 
behind massive traces with a constant one-dimensional density of mass. For the average walking 
path (in the complete multifractal), the geometrical structure built by all traces (weighted with 
the probability of their occurrence) exhibits the same scaling behaviour as the average walk time. 
This geometrical structure will be named "assigned multifractal". Its Hausdorff dimension is 
the walk dimension d w of our interest. 

Just one can question the possibility to use the time scaling factors A found out by trea- 
ting all elementary generators separately. And what happens with a particle returned to the 
starting-point without having reached any end-point, which starts then into another elementary 
generator? Concerning the calculation itself there are no difficulties; but the physical interpre- 
tation has to be discussed carefully. 

As the same walk dimension d w is valid for all kinds of random walk models, we can postulate 
a certain one. This shall be a stationary process taking place in the network after (k + 1) 
construction steps. Further, it shall be an average process, what means - strictly spoken - 
it represents the average behaviour of particles. Then we consider the subcells of the k-ih 
construction step as black boxes characterized by their substructures. After numbering all these 
black boxes, the only well-known thing is that per each time unit Nj particles start from each 
corner point into the j-th black box. The symmetry of this idea causes the arrival of the same 
number Nj of particles from the internal space of the z-tli black box in each of their corner 
points. Introducing the argument of indistinguishability for particles located at the same point, 
we can suppose that the particles starting from any corner point into a related black box are 
the same who have just arrived from this black box. In this picture the particles can never leave 
their cells. This is exactly what we were looking for. 

It remains the proof of renormalizability for the assigned multifractal. Therefore we refer to 
our method for calculating the time scaling factors A. Including the absorbing end-points into 
Q, X^ and X and symmetrizing the problem by carrying out the sum over all possibilities to 
mark the starting-point we get the same total number of passages for all links in an elementary 
generator. Thus, the transition to the next construction step does not disturb our way of 
consideration. 

Now let us remind the multifractal-generator-matrix M. It is known that we can choose any 
partition of the pij satisfying the condition (10). As we would not complicate our calculations, 
we suppose the pij does not depend on the index i. This yields the equipartition pj = l/nj. 
Then the matrix M becomes 



Here we must distinguish exactly the effect and the meaning of the several factors. The gij sort 
the sum and guarantee the correct treatment of each partial sum and nothing more. The bj 
represent the length scaling behaviour and nothing else. And last but not least the nj stand for 
the "mass"-scaling behaviour and for nothing more. 
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For the assigned multifractal we can take over the gij and the bj from the basic multifractal. 
However, the "mass"-scaling factors nj must be replaced by the "mass"-scaling factors of the 
assigned multifractal's elementary generators, the time scaling factors Aj. Eventually we get 

Mx = (m ij ) x = (aij^J ■ (22) 

In this case we have 

Pj = Y_- (23) 



The general form for any partition of the pij is 



Mx=^gij\jj^j = (gij^jPi j b'j) x with J2^ij X jPij = 1 - 



(24) 



The cells of the basic multifractal have no direct meaning for the scaling properties of the assigned 
multifractal, the only direct relation is to the gij. Further, the assigned multifractal's Hausdorff 
dimension d w is larger than the embedding dimension d e . This leads to the imagination of our 
assigned multifractal's proliferation into the (d e + l)-th dimension. 



4 Example 



We choose the two bases 2 and 3, the embedding dimension shall be 2. The numbers of subcells 
in the elementary generators are 3 and 6. Further, we set gij = ^ , Vi, j and pij = pj. Thus we 
get 

2 r-l 3T + I 



M = 

Demanding the largest eigenvalue to be equal 1 leads to the equation: 



39-1 

2 T-1 



61 

3 T + 1 



1 = 



2 r - 



3«" 



-1 3 r+l 
T + 



6« 



Solving the equation yields the Hausdorff dimension Dq = 1.613. . . and the figure 2. Applying 
the /(a)-formalism we get the figure 3. Just we come to the time scaling factors. Therefore we 
need 

" 



Qi 



Further, we are given 
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Figure 3: The global dimension f(a) vs. the local dimension a 
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Figure 4: The Renyi dimension D q for the assigned multifractal 



Now we write down the matrix M\ for the assigned multifractal: 

2 T 3 T 



M\ = — 

2 



59-1 
2 T 



90 
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Demanding the largest eigenvalue to be equal 1 leads to the equation: 



1 

1 = - 

2 



5«-i 



+ 



V 



no 



9-1 



Solving the equation yields the walk dimension d w = 2.323 . . . and the figure 4. There we have 



D. 



d w (b = 3) and = d w (b = 2). Applying the /(a)-formalism we get the figure 5. 



5 Conclusion 

Obviously the calculation of the walk dimension essential for every kind of diffusion works with 
the simplest walk model. The model of the average-walk-time-multifractal has proven to be 
good. For the dimensions there is the relation: Do < d e < d w . Our method can be seen as a 
generalized version of the renomalization method already used to calculate walk dimensions. 
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Figure 5: The global dimension f(a) vs. the local dimension a for the assigned multifractal 

References 

[1] T. C. Halsey, M. H. Jensen, L. P. Kadanoff, I. Procaccia, B. I. Shraiman: Fractal measures 
and their singularities: The characterization of strange sets, Phys. Rev. A 33, 1141-1151 
(1986) 

[2] U. Renner: Dissertation, Universitai Leipzig, 1992, Kap.: Fraktale und Renormierung 

[3] J. Machta: Generalized diffusion coefficient in one-dimensional random walks with static 
disorder, Phys. Rev. B 24(9), 5260-5269, (1981) 

[4] J. W. Haus, K. W. Kehr: Diffusion in regular and disordered lattices, Phys. Rep. 150, Nos. 
5 & 6 (1987), 263-406 



9 



